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Abstract 



We give a definition of asymptotically locally Lifshitz spacetimes, with 
boundary data appropriate for a non-relativistic theory on the boundary. So- 
lutions satisfying these boundary conditions are constructed in an asymptotic 
expansion. We identify the boundary data with sources for dual field theory 
Q-ij operators, and give a prescription for calculating the one-point functions of the 

field theory operators (including the stress tensor) in the presence of arbitrary 
sources. The divergences in these one-point functions can be cancelled by holo- 
pg \ graphic renormalization, adding counterterms which are local functions of the 

£> ' boundary data. 

in 

Tf : 1 Introduction 

The use of gravitational duals to study strongly-coupled field theories [U [2] has pro- 
vided a unique tool which has shed light on a number of important questions, par- 
ticularly concerning thermodynamics and theories at finite density, which are related 
to black holes in the bulk. This work has been extended over the last few years to 
include applications to field theories of interest to condensed matter physics (see [3j H] 
for useful reviews). The application to condensed matter throws up new questions; 
to address these and fully realise the potential of the holographic approach, we need 
to extend our understanding of the bulk gravitational theories and the holographic 
dictionary. An important example of this is the appearance of anisotropic scaling 
symmetries in condensed matter. The low-energy physics near a phase transition 
may be invariant under t — > X z t, x l — > Xx l , where z is referred to as the dynamical 
exponent. We follow |5j in referring to a field theory which has this scaling symmetry 
(and no boost symmetry) as a Lifshitz field theory. Understanding such symmetries 
holographically requires us to go beyond the familiar context of asymptotically anti-de 
Sitter spacetimes on the gravitational side. 
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A holographic duality for these field theories was proposed in [5]. The proposal is 
that the dual of the field theory vacuum has a bulk metric 

df ? = - r ** d1 ? + 1 * ( tf i + L 2 ^, (1) 

where L 2 represents the overall curvature scale, and the spacetime has d + 1 dimen- 
sions, so there are d s — d — 1 spatial dimensions x\j This metric is referred to as a 
Lifshitz geometry. Such a metric can be realised as a solution in a variety of bulk 
gravitational theories with different matter content. In [5], the bulk theory involved 
two p-form fields with a Chern-Simons coupling. A simpler theory with a massive 
vector (which is on-shell equivalent to the previous theory) was introduced in [6J. 
More recently, ([I]) was realised as a solution in string theory in a number of different 
truncations [7J EJ El [lOj [11] (see [12] for earlier attempts). In this paper, we will focus 
on the massive vector theory of [B], which provides the simplest context for studying 
this geometry. The focus is on the differences in structure between geometries of the 
form ([1]) and the AdS case, so the key results will not depend too heavily on the 
specific form of the matter considered. 

In the AdS context, there is a well-developed holographic dictionary, starting from 
the work of p31 E] • This relates asymptotic boundary conditions for bulk fields to the 
sources for dual operators in the field theory description. The asymptotic boundary 
data for the metric is specified by the induced metric on the conformal boundary of 
the spacetime. The results of |T5[ [T6] show crucially that at least in a neighbourhood 
of infinity, there is a bulk solution for any arbitrary boundary metric, so we can 
consider the correspondence for the field theory on an arbitrary background. 

The correlation functions for operators in the CFT are then obtained in the saddle- 
point approximation by considering the variational derivative of the bulk action with 
respect to the sources. When this is applied to the naive bulk action, these correlation 
functions contain divergences. These can be removed by the process of holographic 
renormalization, adding appropriate local boundary counterterms to the bulk action 
to ensure that it defines a good variational principle for the asymptotic boundary 
conditions [T71 ITS] . An elegant approach based on Hamiltonian evolution in the 
radial direction and a functional differential was developed by [HI [20] . 

The aim of the present paper is to develop a general holographic dictionary for 
the Lifshitz case. Starting with [5] , there has been extensive work on the calculation 
of correlation functions from perturbations in the bulk; but so far, this work has 
only considered spacetimes which asymptotically approach ([]]), treating changes in 
the asymptotic boundary data only perturbatively. Here we extend this to consider 
arbitrary boundary data, constructing metrics which approach ([I]) only locally. In the 
AdS case, this generalisation is achieved by considering an arbitrary conformal class 
of metrics on the conformal boundary of the spacetime. In addition to enabling us to 
treat field theories in curved backgrounds, this also provides a more comprehensive 
(and less coordinate-dependent) perspective on the holographic dictionary, allowing 
us to understand its essential features. 



1 We will mostly focus on the case d s = 2. 



In the Lifshitz case, the spacetime ([I]) does not have a conformal boundary, as 
the timelike direction is treated differently from the spatial ones. Thus, while a sim- 
ilar correspondence should exist, we need to develop a new approach to describe it. 
Building on the work of [2TJ (and similar observations in [22]), in section [2] a definition 
of asymptotically locally Lifshitz spacetimes is given, in terms of appropriate bound- 
ary conditions on a set of frame fields describing the bulk geometry. This enables a 
coordinate-invariant treatment of the different scaling of the timelike direction. The 
rescaled boundary values of the frame fields define the boundary data for the ge- 
ometry; they can be thought of as defining a Galilean geometrical structure on the 
spacetime boundary. If we impose a further restriction, the boundary admits a notion 
of absolute time. 

In section [31 the interpretation of the boundary data from the dual field theory 
point of view is discussed. Some general expectations for the behaviour of the stress 
tensor complex in a Lifshitz theory are reviewed. Then following [21] , a prescription 
for calculating this stress tensor complex is given, considering the variation of an 
appropriate bulk action with respect to variations in the boundary data for the frames, 
holding the boundary value of the matter field(s) with tangent space indices fixed. 

In section HI we specialise the discussion to the massive vector model of [B], to 
enable discussion of more dynamical issues. We review esults of [21 J on the solution 
of the linearized equations of motion on the background (pQ) . We discuss some of the 
open problems arising from that work, some of which will be addressed here and some 
of which remain for future work. In particular, we find that a modification of the 
boundary conditions is required for z > 4; we treat this modification in the linearized 
regime, but leave a full exploration of the modified boundary conditions for future 
work. 

In section [51 the solution of the equations of motion in a large-distance expansion 
is considered. We find that asymptotically locally Lifshitz solutions exist for arbitrary 
boundary data for z < 2. For z > 2, there are restrictions on the boundary data, 
setting to zero the sources for the irrelevant operators. There are then asymptotically 
locally Lifshitz solutions for arbitrary sources for the relevant and marginal operators. 
These are key results, providing an analogue of the Fefferman-Graham expansion [15] 
for asymptotically Lifshitz boundary conditions. 

Finally, the problem of holographic renormalization is considered in section EJ We 
obtain counterterms rendering the expectation value of the stress tensor finite for 
arbitrary asymptotically locally Lifshitz spacetimes. The counterterms are obtained 
using an analogue of the approach introduced for AdS in [191 120] , organising the cal- 
culation in eigenvalues of a dilatation operator. This approach efficiently evaluates 
the counterterms; this is illustrated by an explicit calculation of the first few countert- 
erms. For z > 2, it was observed in [21] that there were divergences in components of 
the stress tensor even when the sources were set equal to zero; some of the response 
functions cause divergences in other expectation values. Here we will see that the 
form of the counterterms is uniquely fixed by requiring cancellation of the divergence 
coming from the sources, but that this same term also cancels the divergence from 
the response functions, so that these are "pseudo-non-local" in the terminology of 



Note added: After the appearance of this paper, partial results on holographic 
renormalisation for asymptotically Lifshitz spacetimes were also reported in [2U |25] . 

2 Asymptotically locally Lifshitz spacetimes 

We want to consider spacetimes which asymptotically approach a metric which can 
locally be written in the form ([T]). This condition can be most easily implemented by 
writing the geometry in terms of a set of frame fields. We write the spacetime metric 
as 

ds 2 = VMNe^eW = r, AB e^e^ + (e«) 2 , (2) 

where spacetime frame indices are written M,N = 0,1,... d, while frame indices 
omitting the radial direction are A, B = 0, . . . , d s . We will also sometimes use spatial 
frame indices I, J running over 1, . . . , d s (recall d s = d — 1 is the number of spatial 
dimensions). Similarly spacetime coordinate indices are /U, v = t,Xi,r and boundary 
coordinate indices a,j3 = t,Xi. We partially fix the gauge freedom by choosing a 
Gaussian normal radial coordinate such that e^ r ' = Lr~ 1 dr, and &■ ' have no radial 
components. 

For the metric (£[)), there is a natural set of frame fields where e^ = r z dt, e^ 1 ' = 
rdx 1 . This motivates defining asymptotically locally Lifshitz boundary conditions by 
requiring that the spacetime admit a choice of frames e^ such that as r — > oo, 



ef = r z e^(r,x a ), eg = re^(r,x a ), 



(3) 



where e a (r, x a ), e a (r,x a ) have some finite non-degenerate limits as r — > oo, which 
we will also sometimes call e a (x a ), e a (x a ). That is, we simply replace the coordinate 
frame fields dt, dx % in ([1]) by some arbitrary non-degenerate frame basis e a (x a ). This 
is the most general boundary condition which allows the metric to locally be written 
as (^ in the asymptotic region. 

Note that we have not imposed any specific falloff condition on the subleading 
parts of e a (r, x a ), e a (r, x a ), beyond requiring that they vanish as r — > oo. This is 
because this boundary condition is meant to be purely kinematical, applying to off- 
shell fluctuations considered in the discussion of the action as well as to asymptotically 
locally Lifshitz solutions of any bulk theory. As in the AdS case, in solving the 
equations of motion in section \5\ we will find that the dynamics will dictate specific 
powers of r which appear in the subleading parts of e a (r, x a ), e a (r, x a ) and in the 
matter fields. But the specific fall-offs will depend on the details of the bulk theory 
(again, as in the AdS case), so we do not wish to introduce them in our definition of 
the boundary conditions. 

The boundary frame fields e a (x a ), e a (x a ) define our boundary data; we will 
see in the next section how these provide sources for the the stress tensor complex. 
Geometrically, they define the general background structure appropriate for a non- 
relativistic field theory on the asymptotic boundary of the spacetime. Such data has 
been previously defined in terms of what is referred to as a Galilean metric [261 EZj , 



consisting of a timelike covector (one-form) field, and a degenerate contravariant 
metric orthogonal to the covector. This is equivalent to our frame data, as the 
rescaled timelike frame field er a (r, x a ) provides a distinguished timelike covector on 
the boundary, while the rescaled inverse metric 

$# = lim r 2 g al3 = e iI)a eL (4) 



r— >oo 



V) 



provides a degenerate boundary contravariant metric, which is obviously orthogonal 
to the timelike covector, g a/3 e a = 0. We can also note that the anisotropic scaling 
implies that the bulk freedom to make local Lorentz transformations of the e^ A ' is 
reduced to the freedom to make spatial rotations of the eS ! \ again as we would expect 
for a non-relativistic theory. 

However, in a non-relativistic context it seems natural to require that the boundary 
have a notion of absolute time. That is, the boundary should have a foliation by 
a family of surfaces defining "moments in time". This is not satisfied for general 
boundary data satisfying fl3]): arbitrary e a identifies only a preferred family of curves 
parallel to the distinguished vector field e? y To be able to use e a to identify a 
preferred foliation by surfaces, we need to require that it is irrotational, 

g(°) A rfe (0) = 0. (5) 

This condition can be conveniently expressed in terms of the bulk Ricci rotation co- 
efficients ^bc 1 " 4 , defined by de^ = il BC A e^ A e^ c \ Requiring e^ to be irrotational 
requires the leading term in VL U ° to vanish asymptotically. If we impose this condi- 
tion, we can choose coordinates such that e^ = x(t, x l )dt for some function x(t,x 4 ) 
(at least in an open neighbourhood). The irrotational condition is also a necessary 
and sufficient condition for the spatial frame fields eS 1 ^ to be surface forming. Thus, 
the e^ will describe a non-degenerate ^-dependent curved metric on the surfaces of 
constant t. 

One might then take ([3D and §5§ to define asymptotically locally Lifshitz bound- 
ary conditions. However, we will see in the next section that assuming that e a is 
irrotational corresponds to setting the sources for the boundary energy flux £ l to zero. 
Thus, to be able to calculate correlation functions involving the energy flux, we need 
to consider violations of ([5]) at least perturbatively. In the remainder of the paper, 
we do not in general impose ([5]) (for z > 2 however, satisfying ([3]) turns out to require 
that we set the sources for £ l to zero, which implies we also satisfy ([5])). 

3 Field theory sources 

We now discuss the interpretation of this boundary data in terms of sources for the 
field theory stress tensor. We make the central assumption that as in AdS/CFT, 
asymptotic boundary data for bulk fields should be interpreted as sources for dual 
operators in the dual field theory. In the case of the metric (frame fields), the ap- 
propriate dual operator should be the stress tensor complex, so we first review a few 
elements of the field theory expectations for its components. 



In a non-relativistic theory, the tensor complex should consist of an energy density 
£ and an energy flux £ l , satisfying a conservation equation (in flat space) 

d t £ + diS 1 = 0, (6) 

and a momentum density V{ and a spatial stress tensor II], satisfying the conservation 
equations (in flat space) 

dtVi + djU/^0. (7) 

As in the relativistic case, invariance under the anisotropic scaling symmetry t — > 
X z t, x l — > \x % implies a tracelessness condition zE + IL/ = 0. The key difference from 
the relativistic case is that the momentum density and energy flux in a non-relativistic 
theory are independent quantities. 

In a Lifshitz theory, there are also some differences in the operator dimensions for 
these quantities. These do not seem to have been discussed earlier in the literature, 
so it is worth considering them here. Since £ is the energy density, and the energy is 
dimension z with respect to the anisotropic scaling symmetry, it will have dimension 
z+d s . This is the marginal dimension for a Lifshitz theory, as it matches the scaling of 
the volume element. Thus, energy density is a marginal operator, as in the relativistic 
case. The conservation equation then implies that the energy flux has dimension 
1z + d s — 1. This implies that it is an irrelevant operator for z > 1. Thus, deforming 
the theory by adding sources for the energy flux will change the UV theory, driving 
it away from the fixed point, quite unlike the relativistic case. 

The momentum density V% will have dimension d s + 1, as momentum is dimension 
1 with respect to the anisotropic scaling symmetry. This is relevant for z > 1. Since 
Vi has dimension d s + 1, the conservation equation implies that U J i has dimension 
z + d s , indicating that it is also marginal. This is consistent with the condition 
z£ + LT^ = 0, which requires that H 3 i must have the same dimension as £. 

The fact that Vi is relevant has several interesting consequences. First, it implies 
that we can generate non-trivial renormalization group flows to the IR by adding 
sources for the momentum density. Secondly, the source for this operator will have 
dimension z — 1, so for z > d s + 2, the dimension of the source is greater than the 
dimension of the operator. From the holographic point of view, this suggests that the 
boundary data associated with this source will fall off more quickly at large r than 
the mode giving the source. We will see this explicitly in the linearised analysis in the 
next section. This is reminiscent of the behaviour of scalar fields in the alternative 
quantization of [28J. A similar z-dependent crossover between the source and the 
expectation value was seen for a Maxwell field on a Lifshitz background in [12] . The 
scalar operator V{V % becomes relevant for z > d s + 2 with the fixed source boundary 
conditions, implying that there is then a relevant deformation of the theory which 
preserves the spatial rotation symmetry. The natural endpoint for this flow is a 
Lifshitz theory with fixed Vi boundary conditions. 

How are the components of the stress tensor complex calculated holographically? 
As argued in [21] , the appropriate dictionary is to identify the timelike frame field e^ 
as the boundary data supplying the sources for £, £\ and the other frame fields e^ 
as the boundary data supplying the sources for Vi, II*. More explicitly, assume we 
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have a bulk action S[e^ A \tp] (where ip denotes whatever matter fields we consider), 
which is finite on-shell and provides a good variational principle for our boundary 
conditions, so SS = for variations satisfying Se^ = Oo Then writing 

5S= [e(T a B 5e^ + 0^) } (8) 



where e is the coordinate-invariant volume density given by e = e^ Ae^ A . . . Ae^ da \ we 
identify T a with the vacuum expectation value (vev) of the energy density £ and the 
energy flux £ *, and T°j with the vev of the momentum density Vi and the stress tensor 
IP . In accordance with the usual holographic dictionary, we also identify O^ with 
the vevs of operators dual to the matter fields. Note that to have the Lifshitz solution 
([I]), the bulk theory will need to include vector or tensor fields which are non-zero 
in the Lifshitz geometry, which dynamically select the timelike direction as different 
from the spacelike ones. Following [29] [2T] the ifi are then understood to include the 
tangent space components of these vector or tensor fields. That is, the variation of 
the frame fields in (JSJ) is with the tangent space components of these quantities fixed, 
rather than the spacetime components. This is an essential ingredient to ensure that 
T l and T t I are different, so that they can be identified with the distinct components 
of the stress tensor complex [30]. It can also be understood physically by arguing 
that we take the variation with the tangent space indices fixed so that the alignment 
between the distinguished frame field e^ -* and the matter fields that are responsible 
for singling it out is maintained as we do this variation. 

One simple check of this prescription is that these quantities are conserved as re- 
quired from the field theory point of view. If we assume that the action S is invariant 
under diffeomorphisms of the boundary coordinates, then as shown in the asymp- 
totically AdS context in [29], this diffeomorphism symmetry implies a conservation 
equation, 

V a T a p - 0^7 ^ = 0, (9) 

so T a p is conserved up to the presence of sources. This is the analogue of (jH]), © for 
general boundary sources. 

Note that the sources for the energy flux £ % are e\ , the dx % components of the 
timelike frame field. Thus imposing the irrotational condition, which allows us to 
set e\ =0 by choice of coordinate system, is setting these sources to zero (up 
to diffeomorphism). Since £ l is an irrelevant operator in the field theory, it is not 
surprising that allowing arbitrary sources for it leads to a qualitative change in the 
UV behaviour, violating what we may want to call asymptotically Lifshitz boundary 
conditions. The precise form of this change is however a little unusual. At least for 
small z, it is not that some field is growing too quickly (as we will see again in the 
next section), but rather there is an obstruction to defining an absolute time on the 
boundary. 



2 This is a non-trivial assumption, but section [6] shows that we can provide such an action by 
holographic renormalization. 



4 Massive vector theory 

Our discussion so far has been very general; we will now turn to the specific example of 
the massive vector theory of [6], in four bulk spacetime dimensions (so d s = 2), so that 
we can consider more dynamical issues. In this section, we will give the definition of 
this theory and review the previous work of [21] on the linearized equations of motion 
on the background (1T1). We will also comment on issues to do with the boundary 
conditions in light of this linearised analysis. In the next section we construct solutions 
for arbitrary boundary data in an asymptotic analysis. While we focus on the massive 
vector model as a specific example, the extension of this analysis to for example the 
string theory truncations of [3 El QUI E] is in principle straightforward. 
The bulk spacetime action for the massive vector theory is 

i r «, , ~. i„ ™„ i o i 



S = ——— I <Pxy/=j(R-2A-^-F ^F^~m 2 A ^)-—- I d^V-hK, (10) 

107rG4 J 4 2 07TG4 J 

where we have included the Gibbons- Hawking surface term, so that this is a well- 
behaved action principle for manifolds with boundary with Dirichlet boundary con- 
ditions. The equations of motion for this theory are 

R„u = A<V + \f^xF u x - ^F Xp F x fg^ + ^m 2 A p A v (11) 

and 

V M F^ = m 2 A v . (12) 

To have the solution ([TJ, we choose A = — ^(z 2 + z + 4) and m 2 L 2 = 1z. Then the 
theory has a solution with metric (CQ) and 



A = ar z dt = J^ — l X z dt. (13) 

4.1 Linearized analysis 

The linearized analysis in |2Tj found a general solution of the linearized equations 
of motion around the pure Lifshitz background §B)U From this linearized analysis, 
one can identify the modes corresponding to an infinitesimal change in the sources 
given by the boundary values of e^ A \ and the modes corresponding to the part of the 
solution of the equations of motion which is not locally determined in terms of these 
boundary values. The essential features can be seen from the analysis of linearized 
perturbations which are constant along the boundary directions. Written in terms of 
frame fields, the results of [21J are that 

e (o) = r z^ + rj dt + rvj . dx i i e W = r z v ^ dt + r ($i + lk%)dx j , (14) 

2 J 2 J 



3 A similar analysis also appeared in [31]. Recent extensions appeared in [32, 24 



where 

/ = c 4 + ^ Cl r^+ 2 ) + 2 (5^-2-l) i ( , +2+fe) _ (5s-2 + &) i ( , +2 _ fe) 
2 + 2 (2 + 2 + /3 2 ) (2 + 2-&) 

(15) 

vii(r) = c li r( 2 - 1 )+c 2i r- 3 + c 3j r-^ +1 ), (16) 

^(r) = C4 . r (i-) + lf!^ C2 . r -3 + ^^ C3 . r -(2 Z +i) 
v ; z(z-4) (2 + 2) 



and 



with 



fcy — dijK + fcj • , (.1' J 



(^ + 2) (z + 2 + f3 z ) (2 + 2-/3,) 

(18) 
and 

C T = C = *«n + ter-W, k T J y = k T J = t ol + t o2 r-^ +2 \ (19) 

The massive vector field was taken to be of the form Aa = S a Aq by choice of frame, 
choosing the frame vector eus\ to be aligned with the part of the vector field parallel 
to the boundary, and A = a + j, with 

= _(£±l) Cir -(,+2) _ (i±iL 2r -i(*+2+A) + Ii±l) C3r -i(.+2-/3,)_ (2Q) 

(2-1) (2-1) (2-1) 

Here /3J = 9z 2 — 20z + 20 = (z + 2) 2 + 8(z — l)(z — 2). Note these results are for generic 
values of z. For specific even integer values, there will be logarithmic solutions; see 
[2Tj for details. Also, here we have redefined vu, v 2 i compared to [21J to highlight the 
field theory interpretation of the bulk modes. We adjusted the explicit powers of r in 
front of Vu, v 2i in e^ so that all of /, j, Vu, v 2i , ky are unchanged under the dilatation 
isometry r — > A _1 r, t — > X z t, x — > Xx of the background (JT]). This ensures that the 
scaling dimension of the different modes in /, j, Vx%, v 2 i, kij can simply be read off from 
the powers of r associated with them. 

This is a solution of the linearized equations of motion for arbitrary constant 
values of the coefficients. If we promote the coefficients to functions of the bound- 
ary coordinates, these are still solutions to leading order in an expansion in r, but 
there will be further subleading terms determined in terms of the derivatives of the 
coefficients, as described in detail in [2T] . 

The physical significance of these modes can be readily identified: 

• The modes C4, cu are changes of lim^oo &°\ Accordingly, these are interpreted 
as sources for £, £ l . 

• The modes c 4i , c 5 , tax, t x are changes of lim^oo eS 1 ' . Accordingly, these are 
interpreted as sources for Vi, H Z j. 

9 



The modes c\, t^2, £ D 2 have scaling dimension z + 2. They were shown to give 
finite contributions to the expectation values of S, IP in [2T] . 

The modes C21 have scaling dimension 3 and were shown to give finite contribu- 
tions to the expectation values of Vi in [21] . 

The modes c^ have scaling dimension 2z + 1 and were shown to give finite 
contributions to the expectation values of expectation values of S l in 



The remaining modes C3, C2 can be interpreted as the source and expectation 
value of an operator O^ associated with the massive vector field. From their 
scaling dimensions, we can read off the dimension of this operator, A^, — Mz + 
2 + j3 z ). Note that this operator will be relevant for z < 2, and irrelevant 
for z > 2. As one might expect for an irrelevant mode, the source C3 makes 
contributions to e^°\ e^ which violate the boundary conditions when z > 2. 
The mode c<i was shown to give finite contributions to the expectation values 



of Oj, in [21]. 

We see that the solutions of the linearized equations of motion divide into source 
terms and the corresponding response functions, as expected. Note that because 
of the conservation equations and the conformal constraint relating energy density 
to spatial stress, the number of independent response functions and the number of 
sources does not match up. This could be resolved by working with gauge-independent 
combinations of the sourcesjj 

One might have expected that we should have a vector operator corresponding to 
the field A^, but this is not correct: because we can choose Aj = by a choice of 
frame, there is only one piece of boundary data associated with this field, A , so there 
is a single scalar operator dual to changes in this boundary data. Physically, what 
happened is the spatial vector worth of additional boundary data in A^ was absorbed 
by the frame fields, providing the extra information needed to have independent 
sources for £ l and V%. 

There are a few interesting open issues already at the level of the linearised anal- 
ysis. In [21], the prescription of section [3] was shown to give finite expectation values 
for all components of the stress tensor complex, once we added appropriate coun- 
terterms to the action. Surprisingly, there were divergences in £ l even when we set 
the modes associated with changes in the boundary data to zero, which required 
some derivative counterterms. In section [6j we will see that the counterterms for 
asymptotically locally Pifshitz boundary conditions are uniquely fixed by cancelling 
the divergences arising from the boundary data. These counterterms should then also 
cancel the divergences involving the response functions; we comment on this issue in 
the discussion. 

There was also an issue with a divergence in the expectation value of O^ for z > 2. 
We will not investigate this further here, as this operator is irrelevant for z > 2, so 
considering sources for it would take us outside the asymptotically locally Pifshitz 



4 An error in the first version of this paper has been corrected here in light of the discussion in 
[24]. 
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boundary conditions we are considering in this paper. The analysis could simply be 
extended by considering this source perturbatively, but we leave this for future work. 
However, we should note that there is an additional issue here which does not appear 
for example for irrelevant scalars in AdS: because of the coupled structure of the 
equations, the source mode C3 appears linearly in the frame fields as well as in the 
vector field. Thus, we cannot think of varying C3 as simply a variation of A , even 
perturbatively; it necessarily involves a variation of the frame fields as well. 

The most interesting open issue, however, is the crossover between the source 
mode c 4 j and the expectation value c 2 j in v^%- When z > 4, the linearized mode 
with coefficient C21 falls off slower at large r than the c^ mode. As a result, turning 
on non-zero c<n will violate our asymptotically locally Lifshitz boundary conditions. 
That is, imposing asymptotically locally Lifshitz boundary conditions for z > 4 gives 
us a smaller space of solutions than expected in the asymptotic regime; this will 
generically lead to problems satisfying the regularity conditions in the interior of the 
spacetime for arbitrary sources. We are therefore motivated to find different boundary 
conditions in this regime. At the linearized level, it is easy to see what we should do: 
considering the mixed boundary condition 

lim r- x {e { p - {z - 3)eJ 7) ) = (4 - z)v\x a ), (21) 

r— >oo 

v\x a ) would provide a source for the momentum density "P% while leaving the mode 
with e\ ~ r z ~ 3 free. Or we could continue to impose a Dirichlet boundary condition, 
but for z > 4 relax this to require that at large r, 

e { p ~ r z -' i w\x a ). (22) 

The mode with e t ~ r is then free. This should correspond to a Legendre-transformed 
theory where we regard V{ as a source current, coupled to a vector operator of di- 
mension z — 1. The value of W7 7 (x a ) then fixes the value of Vi, thought of as the 
source. The Dirichlet boundary condition fixes the leading term in e t , so it is the 
analogue of the usual quantization for a scalar, while the mixed boundary condition 
fixing the source for V 1 is the analogue of the alternative quantization. As noted in 
the previous section, for z > 4 the operator ViV 1 generates a flow from the fixed 
source boundary condition in the UV, which we expect to run to fixed V % boundary 
conditions in the IR. Hence the latter are more generic in this range. Extending this 
analysis beyond the linearised level is an important open question, but we will leave it 
for future work and focus on the asymptotically locally Lifshitz boundary conditions 
(J3J) in this paper. 

5 Asymptotic expansion 

In this section, we want to go beyond the linearised analysis by showing that solutions 
of the bulk equations of motion exist for arbitrary boundary data. To do so, we will 
solve the equations of motion in an asymptotic expansion: that is, we work at large 
r, and solve the equations in an expansion in powers of r. To treat this large r 

11 



expansion, it is convenient to rewrite the equations of motion in a radial Hamiltonian 
framework. If we work in Gaussian normal coordinates, the canonical coordinates on 
a surface r = r are the induced metric h a/3 and the gauge field A a . The momenta 
are 7r a/3 = K a p — h a pK and TT a = F ra , where K a $ is the extrinsic curvature of the 
surface Jj 

To make contact with the previous section, we may note that the momenta 7r Q/ g, 
Tr a can be related to the expectation values of the operators in the dual field theory 
by observing that on-shell, the variation of the action (fTUj) is 

5S= [cPxV^iir^Shais + irSAa)- (23) 

When we do frame variations Se a holding Aa fixed, the variation of the action 
on-shell is then 

5S = I d 3 xV^h((27T a p + 7T a A p )ef B) Se^ + tt a 5A a ), (24) 

where ir A = Ti a e a • Thus, we can define an object related to the boundary stress 
tensor complex considered in the previous section: 

T^ = (2n a , + n a A^ {B) . (25) 

Note that the variations in ( 124"|) are with respect to bulk fields, whereas the variations 
in (jSJ) are with respect to the boundary data. Thus this bulk object will differ from 
the boundary stress tensor by some powers of r. We also call this the "bare" stress 
tensor because using the naive Dirichlet action ( TlOl) gives a divergent stress tensor 
which will require renormalization. We address this in the next section. The point 
at this stage is just to note the relation between the momenta and the stress tensor 
complex. Unlike in the AdS case, we cannot simply trade the stress tensor for n a p; 
the T aB also include information from 7r a . 

5.1 Gauss-Codazzi equations 

To solve the equations of motion in an expansion in powers of r, it is useful to 
decompose the equations in the radial direction, using the Gauss-Codazzi equations 
to rewrite them as a pair of dynamical equations, 

K a(3 + KK a/3 - 2K ai K~<p =R a p - Ah a p - -F ai Fj + -h a(} F la F J ' 7 - -n^p 

+ -hapn^Tr 1 - -m 2 A a Ap, (26) 

n a + Kn a + VpF 130 = m 2 A a , (27) 



5 Canonical momenta would include a factor of y—h, but this convention will simplify the relation 
to the vevs. 
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V a K% - VpK a a = 


: -F Pa Ti a + -m 2 ApA 


K 2 


- K afi K^ 


= R-2A + ~7r a 7r Q - 


~ \Faf>F* + \m 2 A 2 T 


and 









and a set of constraints, 

1 1 _ 

(28) 

\m 2 A a A a , (29) 

V a vr a = -m 2 A r . (30) 

Here ' denotes rd r , and the Ricci tensor R a p and covariant derivatives V a are with 
respect to the induced metric h a p on a surface of constant r. Using fT25|) and (j30|) . 
we can see that the constraint equation f )28|) encodes the conservation equation (J9J. 
This is as expected; this is the diffeomorphism constraint in the radial Hamiltonian 
framework, associated with coordinate transformations on the surfaces of constant r. 

With our boundary conditions, the constraints (|28|) and (!29l are automatically 
satisfied when the dynamical equations are. This is because in general, the dynamical 
equations imply the r-derivative of the constraints, so the only non-trivial additional 
information in the constraints is the r-independent term. With asymptotically locally 
Lifshitz boundary conditions, the r-independent term reduces to the constraint in 
a Lifshitz spacetime, which is satisfied because this is a solution of the equations 
of motion. Thus, we need only consider the dynamical equations f )26|) and (T27|) . 
substituting f )30|) for A r . 

Since the asymptotic boundary condition is written in terms of frame fields, it will 
be convenient to also rewrite these equations in terms of frame fields e^ A \ As in the 
linearized analysis, we partially fix the choice of frame by setting Aj = 0, so that the 
distinguished timelike frame field is always aligned with the boundary component of 
the massive vector field. Note that the massive vector field will still have a non-zero 
radial component, but this is determined algebraically by ( 130]) . The relevant data is 
then the frame e a and the matter field Aq. 

For convenience, we define a "frame extrinsic curvature" 

K A B = ef B / a A \ (31) 

Note that this is not a symmetric object, and the usual extrinsic curvature is K a p = 

e a ei K(ab)> where the round brackets denote symmetrisation. The equations of 
motion written with frame indices are then 

K {AB) + KK (AB) + \{K CA K C B - K AC K B C ) =R AB - Arj AB - -F AC F B C + IvabF CD F cd 

1 1 c 

- -jK A k b + -Vab-kck , (32) 

v A + Kn A - K A B n B + V B F BA = m 2 A A , (33) 

and the constraints 



\k a <k a = R-2A- \f ab F ab + \m 2 A 2 r - I, 



K 2 - K {AB) K AB - -ti a ix a = R-2A- -F AB F AB + -m 2 A 2 r - -m 2 : A A A A , (34) 
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V A K (AB) - V B KJ = l -F BA ir A + l -m 2 A B A r) (35) 

V a tt a = -m 2 A r . (36) 

Here Fab = e ?A) e (B)Fa/3, and V^ = e^Va, where the covariant derivative V Q is a 
total covariant derivative (covariant with respect to both local Lorentz transforma- 
tions and diffeomorphisms) , determined by requiring V^el = 0. The connection 
Wage defined by 

V a V B = d a V B - u aB c V c (37) 

is then related to the Ricci rotation coefficients by [33] 

LO ABC = —&ABC + &ACB + &BCA, (38) 

where as a reminder we define the Ricci rotation coefficients by de^ = £lAB° e 'Ae*^ '. 
This implies uj, ab ^ = —VL A b C i an d w^d = ^cd- Note that VL AB C * s antisymmetric 
in its first two indices, but the orthonormality of the frame fields implies that ojabc 
is antisymmetric in its last two indices. 

5.2 Asymptotic analysis 

We want to study the equations of motion in the asymptotic regime by making an 
expansion in powers of r. The aim in this section is to analyse the structure of 
the equations of motion and demonstrate that a solution exists. We will find that 
asymptotically locally Lifshitz solutions exist for arbitrary values of the sources for 
z < 2, and if we set some of the sources to zero solutions continue to exist for z > 2. 
We first consider z < 2 in detail, and then comment on the differences for z > 2. 

We recall that e^ ^ = r z e^ and e^ = re^\ and write A = a + ifj. The leading 
term in e^ A ' is independent of r, and gives the boundary data dual to the stress tensor. 
From the linearized analysis, we see that ip = r~ A -i/j(r, x a ), where A_ = z + 2 — A^ = 
\{z + 2 — (3 Z ). The asymptotic value lim^oo if) corresponds to the source for the dual 
operator O^. For z < 2, this is a relevant operator, so A_ is positive. We want to 
show that having chosen these leading terms, we can iteratively find a solution of 
( 1261) . ( |27|) in an expansion in powers of r by adding subleading terms to &• > and ijj. 

It might appear that we have fewer equations than unknowns, since we have nine 
components e a plus one ip, and only nine equations. This is because there is a 
residual freedom to make a local rotation in the choice of the spatial frame vectors 
e^\ giving us one pure gauge modejj It is convenient to fix this gauge symmetry by 
choosing the spatial frames so that Kjj is symmetric. 

Although e^ and ip are the fundamental quantities we are solving for, in analysing 
( 1261) and (1271) it is more convenient to work in terms of Kab and ii A . These are given 



6 In general dimension, there are d 2 + 1 fields, but the spatial rotations make ^d s (d s — 1) = 
\{d— l)(d — 2) of them pure gauge, leaving \d{d + 1) + d physical fields, matching the number of 
equations. 
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by 

^oo = -^ + e ( o) a e" 0) , (39) 

Kqi =r z ' 1 e { o) a ef I)j 
Kio =r 1 - z e {I)a ef 0)l 

which implies K = z + 2 + e(A)a& a j and 

7T = i> - Koo(a + ^) - <9 A r , (40) 



ttj = - K 0I (a + ip) — djA r . 



Working in an expansion in powers of r, we can solve (136]) at each order to determine 
A r . Because of the derivative, this will be determined in terms of contributions at 
earlier orders. 

There is an explicit positive power in Kqi', to have a solution where all quantities 
involve only negative powers, the dependence of e(o)c*e?U must be sufficiently sup- 
pressed so that the expansion of Kqj contains only negative powers. We will see that 
we can self-consistently find such a solution. 

The curvature and field strength terms will also involve explicit powers. The key 
point is to see that these all involve negative powers, so that the sources produced 
by considering arbitrary boundary data can be cancelled by adding subleading terms, 
and the corrections produced by the curvature of the subleading terms are further 
subleading. This turns out to be true for the Ricci rotation coefficients: 

V=r-Ve (0) We[o)<J)> ( 41 ) 

\ljj =r [de )a/3 e (7) e (j)> 

"(XT ~~ ' \ at - /«/3 e (0) e (7)' 

\ljj =r (de )ap£(i}&(jy 



Explicit derivatives will also come with similar factors: d = e? Q ,d a = r z e? sd a and 
<9 7 = ef I} d a = r -1 e^a a . 

We might have expected some positive powers to appear, since we are considering 
turning on sources for the irrelevant operator S\ and the breakdown of the UV be- 
haviour associated with irrelevant operators is normally signalled in the bulk through 
the appearance of growing terms in these equations (see [23J for a recent discussion 
in the AdS case). The term that comes closest to this expectation is fi/j°, but it 
has an extra factor of r _1 because it involves the derivatives of the source for this 
operator, rather than the source itself. A constant value of this source corresponds 
to a diffeomorphism, so only its derivatives will enter into the equations of motion, 
and thus for z < 2 it will not lead to positive powers in the expansion despite the 
operator being irrelevant. There is still a change in the UV behaviour when we have 
non-zero sources for S l , but it is signalled only by the violation of (jSJ), implying that 
the boundary does not have a notion of absolute time. 
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Using (1351) . we can express the curvature terms in the field strength and the Ricci 
tensor in terms of the Ricci rotation coefficients. For the field strength, 

F AB = d A A B - d B A A - uj ab c A c + oo BA c A c = d A A B - d B A A + 2Q AB C A C . (42) 

Hence Fqi = —diAo + 2Q qi Aq will contain an explicit r" 1 , and Fu = 2Qjj°Ao will 
contain an explicit r z ~ 2 . The covariant derivative term appearing in (133)) is 

V A F AB = d 2 A B - d A d B A A + 2d A (n AB c A c ) - 2n c DC F DB - n AD B F AD , (43) 

so \7 a F A q will have terms containing r~ 2 and a term quadratic in fijj containing 
r 2z ~ 4 , while V A F A j will have terms containing r~ ( ' 1+2: - ) and a term linear in Vtjj 
containing r z ~ 3 . The Ricci tensor is 

Rab = - Qa^cb ~ dcu AB c + u CD au DC b + u ABD u)(? D , (44) 

= — 2d A Q CB + dc{£l AB + Q B A + Q A B ) — £1 C da&> b + ^cad^ b 

— 2Q C {Vt ADB + Qbda + ^adb)- 

As a result, i?oo an d Rij will contain terms with r~ 2z , r~ 2 , and a term quadratic in 
fljj containing r 22-4 , while Rot will contain terms with r~( 1+z ) and a term linear in 
Jlfj containing r^ z ~ 3 \ 

Thus, all the explicit powers are negative, and we can construct a solution by 
making an expansion in negative powers. Let us write 

K AB = J2r- A K { A A l ^ = £r-V A) . (45) 

A A 

Assuming the sum only involves negative powers, that is A > 0, then the terms in 
(|26|) . (|2"Tj) of order r _A ° will involve K AB \ tp^ A ' for A < A , and curvature terms. 
The 00 and IJ components of (|26j) and the component of (l2"Tj) can be solved to 
determine the A = A components of -Koo> Kjj and ip, while the 01 components of 
( 1261) and the I components of ( 1271) determine the A = Ao components of Kqi and 
Kjq. The sources are given by the quadratic terms involving the extrinsic curvature 
terms and tt a at earlier orders, and the curvature terms which involve the frame fields 
and ip at earlier orders. 

To determine which powers occur, we work iteratively starting from the leading 
source terms coming from derivatives of the boundary data. The curvature of the 
r-independent terms in e a will give sources in the 00 and IJ components of (1521 
and the component of (127)) . These sources go like r~ 2z , r~ 2 , and r 2z_4 . Similarly 
substituting the curvature of the r-independent background values will produce terms 
in the 01 components of (!32j) and the / components of (!33j) which go like r~( 1+2 ) and 
r z ~ 3 . To cancel all of these terms, we need subleading terms in K AB , 4> with these 
same falloffs. Thus, we would need subleading terms in the frame fields such that 

e (o)« e (o) ~ r ,r ,r , (4b) 

ki)^ a {J) ~r 2z -\r- 2 ,r- 2z , 
e(o)aef 7) ~r~ 2 ,r~ 22 , 
kifrtfk ~r 2z -\r~ 2 ,r- 2z , 
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Note that because of the explicit r z 1 in K 0I , the required behaviour for e( ) a e?U does 
not involve the leading r 2z ~ 4 falloff that appears in all the other terms. That is, the 
presence of the explicit positive power in K j has the effect that the behaviour of the 
subleading terms in e(o) Q obtained by solving these equations will produce a faster 
falloff for the term contributing to Kqi. 

If we include the effect of the leading term in tjj, this will give contributions to the 
equations with powers of r~ A ~. This can appear on its own in the 00 component of 
( 132]) and in the component of ( 133]) . and in combination with derivative terms in all 
the equations. Since Aa appears quadratically in (13"2"]) . there will also be terms going 
like r~ 2A ~ . To cancel these contributions, we will need to add additional subleading 
contributions containing powers r~ A ~. These are negative powers for z < 2. 

All of these subleading terms will then in turn contribute to the curvature terms 
and to the quadratic terms in ( 126]) . ( 127"]) requiring further subleading terms to cancel 
these contributions. We will thus need an infinite series of contributions to e^ A \ ip, 
with powers 

A = 2m + 2nz + (2z - 4)p + A_g, m,n,p, q E 0, 1, 2, . . . (47) 

We have concentrated above on subleading terms in the large r expansion which 
are locally determined by the boundary data. However, there will also be independent 
terms in the expansion which are not locally determined by the boundary data, cor- 
responding to the second set of linearised solutions in section 14.11 These are referred 
to as the response functions, as they determine the finite part of the expectation 
values, as can be seen from the linearised analysis of [21J. These appear first at spe- 
cific characteristic powers of r corresponding to the dual operator dimensions. The 
leading-order terms will be given precisely by the linearised solution in the previous 
section, as the x a dependence of the boundary data does not effect the radial equation 
for the leading-order part of these modes. 

In the asymptotic expansion, we can see the leading r-dependence of these addi- 
tional contributions by finding the value of A where the equations at order A have 
a non-trivial solution with no source contribution. This is most easily illustrated by 
considering the trace- free part of the IJ components in (126]) . for which the leading 
contribution to the equations at e~ Ar is simply —AK}j + (z + 2)K}j , so that in 
the absence of sources there is a solution at A = z + 2, corresponding to the modes 
td2 and t 2 in section 14.11 This works similarly for the other modes, but the coupling 
between the equations makes the analysis more complicated. For generic z, these 
response functions appear at a value of A where there is no source contribution. For 
some z, the A determined in this way is one of the values ([47]) . and then there will be 
a logarithmic solution; we will ignore these cases for simplicity, but it should be easy 
to include them. The curvature terms arising from the leading part of the response 
functions will source another infinite series of subleading corrections. 

Thus in summary, for z < 2, there are solutions of the bulk equations of motion 
for arbitrary boundary data, which can be constructed order by order in an expan- 
sion in powers of r. There is an infinite series of terms with powers ( j4"Tj) which are 
determined by the boundary data. The solution in this asymptotic expansion is not 
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unique, as there are additional response functions in the general solution which are 
not determined locally in terms of the boundary data. One expects these to be fixed 
in terms of the boundary data by imposing appropriate regularity or inner boundary 
conditions in the interior of the spacetime. In our analysis of the counterterms in the 
next section we will assume this is true, so that the bulk solution is uniquely deter- 
mined by the boundary data, but we will not explore it in detail, as we are restricting 
to an asymptotic analysis. 

5.2.1 Expansion for z > 2 

In the above analysis, we assumed that z < 2. If z > 2, then A_ < 0, and to have an 
asymptotically locally Lifshitz solution, we need to set the source term lim r _ s . 0O ifi to 
zero. This is a familiar story from the AdS context; adding a source for an irrelevant 
operator will change the UV behaviour of the theory, which is reflected in the bulk 
dual by a violation of the usual asymptotic boundary conditions. 

Additionally, the explicit power appearing in fi 7 j° is now positive, so in order to 
have a solution with all negative powers, we have to restrict the leading contribution 
to e a such that it is irrotational. That is, we now need to set the source for the 
energy flux E l to zero to have asymptotically locally Lifshitz solutions satisfying fl3]). 

With this restriction, the leading contributions in the 00 and I J components of 
(I52"|) and the component of (j2"7|) from the curvature of the boundary data will go 
like r~ 2z and r~ 2 , and the leading contributions in the 0/ components of (|32|) and the 
I components of (I33j) from curvature of the boundary data will go like r~^ l+z \ To 
cancel these, we need subleading terms in Kab, ip with these same falloffs. Thus, we 
would need subleading terms in the frame fields such that 

%) )a ef 0) ~r- 2 ,r^, (48) 

e(i) a efj) ~r~ 2 ,r- 2z , 
e(o)aefo ~ r~ 2z , 
e(7>e (0) ~ r , r 

Again, we see a faster falloff in e(o) a e?U. As before, this faster falloff ensures that K 0I 
will involve the same powers as Kj , despite the positive power in the expression for 
Kqi- It also implies that the first non-zero term in fljj comes from the terms at 
A = 2z, and hence goes like r~ 4 . 

Thus, for z > 2, once we fix the sources for O^ and £ l to be zero, the powers 
appearing in the expansion of (I52"j) and (1271) will start with r~ 2 and r~ 2z , and for 
arbitrary values of the other boundary data, there will be a solution in a series of 
negative powers as before, with powers 

A = 2m + 2nz, m,neN, (49) 

together with the response function contributions from section 14.11 and their sublead- 
ing corrections. 



Note however that as noted in section I4.1[ if we impose asymptotically locally 
Lifshitz boundary conditions for z > 4, the size of the space of response functions is 
reduced, because the mode c 2 i has to be set to zero to satisfy the boundary conditions. 
We therefore expect that we will not generically have enough freedom to satisfy the 
inner boundary condition for arbitrary boundary data. That is, for z > 4, while we 
will have asymptotically locally Lifshitz solutions in the asymptotic regime, to have 
an appropriate space of regular solutions for arbitrary boundary data we will need to 
consider a different boundary condition as outlined in section 14.11 

6 Holographic renormalization 

Having explored the space of solutions for our boundary conditions, we should now 
construct an appropriate action principle for arbitrary boundary data to complete 
the holographic dictionary. To be able to apply the prescription given in section [31 
we need an action which is finite on-shell and has vanishing variation on-shell under 
arbitrary variations satisfying the asymptotic boundary conditions. That is, SS = 
on-shell if the field variation does not change the boundary data. We will construct 
such an action by adding appropriate local boundary counterterms to cancel the 
divergences in the action (TTOj) . 

We will give an efficient algorithm for obtaining the required counterterms, based 
on the approach of [I2J [20] • This involves expanding in eigenvalues of an appropriate 
bulk dilatation operator, rather than in powers of r as we did in the previous section. 
The first step is to take a Hamilton- Jacobi style approach to determining the on-shell 
action. Assuming that we impose some appropriate boundary or regularity condition 
in the interior of the spacetime, the on-shell solution of the equations of motion will 
be uniquely determined in terms of the asymptotic boundary data, so the on-shell 
action is a function of the boundary data, which we can write as a boundary term, 

S= [ d 3 xV^h\{e (A) ,ip). (50) 



Note that we actually write the action as a function of the values of the bulk fields 
evaluated on the boundary, rather than in terms of the boundary data. For asymptot- 
ically locally Lifshitz boundary conditions, the leading asymptotic values of the fields 
are determined by the boundary data, so this contains equivalent information, but it 
is better to work in terms of the boundary values of bulk fields, since counterterms 
added to the action need to be written in terms of the latter. 

The variation of the on-shell action is given by f l23|) . and we can use this to replace 
the dynamical ODEs fl32l) . fl33|) by functional differential equations, writing 

*"* = */ S, t^ = _L *5. (51) 

yJ-hoh a p yJ-hoA a 

These relations give the canonical momenta as functions of the boundary data. As- 
suming that the on-shell action is boundary diffeomorphism invariant will imply that 
these quantities automatically satisfy the diffeomorphism constraint (1281) . and if we 
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substitute them into the Hamiltonian constraint (|29|) . this becomes a quadratic equa- 
tion in the functional derivatives of A, which can be solved to determine A in an 
appropriate expansion. The divergences we need to cancel are determined by apply- 
ing this procedure to the naive action (ITU]) , calculating its on-shell value as a function 
of the boundary data. Since we are working with divergent quantities, we will need 
to introduce a regulator. We therefore assume that the boundary integral in ( 150]) is 
evaluated on a surface r = r , and the value of the on-shell action is then a functional 
of the values of the bulk fields evaluated on this surface; these Dirichlet data deter- 
mine the value of the action for the bulk solution. Also note that the value of the 
on-shell action is in general a non-local function of the boundary data, as it involves 
the response functions, which are only determined once we impose regularity condi- 
tions in the interior of the spacetime. An essential point is to see that the divergent 
terms in the naive on-shell action f llOp are local functions of the boundary data, and 
do not involve the response functions, so they can be cancelled by simply subtracting 
the same function as a boundary counterterm in our definition of the action. 

For the Lifshitz case we need to introduce slightly different variables, defining the 
functional derivatives 




(52) 



and 



*V> = -^f4t S = -^^-TT^TT [ d 3 xV^h\. (53) 

This change of variables is important because our boundary conditions give definite 
falloff behaviour for the frame fields, not the boundary metric. From the field theory 
point of view, (152"]) and ( 153"]) are convenient variables because they will determine the 
vevs of the corresponding boundary operators (we explain the relation to the vevs 
in a little more detail below). We will choose a gauge such that Tjj is symmetric, 
so that the number of quantities in (1521) and (|53|) matches the number of equations 
we are replacing. From the relation (1251) . we see that we can obtain the canonical 
momenta from these new variables, by writing 7r° = n^, using 

ttjAo = T I0 - T 0I (54) 

to determine 717, and recovering 11 ab from 

kab = 2^ Tab ~ ^aAb)- (55) 

Thus, the canonical momenta are determined in terms of functional derivatives of 
the action, and the constraint ( 134"|) can be solved in an appropriate expansion to 
determine the value of the on-shell action density A. 

We want to introduce a functional derivative which agrees at leading order with 
the radial derivative, which will be used to organise the expansion of the on-shell 
action. Assuming our asymptotically locally Lifshitz boundary conditions (J3j), we 
have that e a at leading order goes like r z and e a at leading order goes like r, while 
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from the linearized solution, at leading order if) goes like r A ~ . Thus, if we introduce 
the dilatation operator 

fo = / A f>^ + e<"^ - (z + 2 - A, W >£) , (56) 

then acting on any function of e^, ip, this will agree with the radial derivative to 
leading order at large r, So ~ rc? r J3 We can also see that <5d can be viewed as a bulk 
analogue of the field theory dilatation operator; the field theory operator would have 
the same form as ( | 56 [ ). but with the bulk fields replaced by boundary data. From 
(152|) and (1551) . we can see that if we consider the dilatation operator acting on the 
on-shell action, we can replace the functional derivatives on the RHS acting on the 
action by T A B , tc^. As a result, we have a relationship 

(z + 2 - 5 D )X = zT° + T I I -(z + 2- A^vr^. (57) 

We now expand all these quantities in eigenfunctions of the dilatation operator 
5d, writing 

A = ^A( A ), S D X^ = AX^\ (58) 

A>0 

and similarly for T A Bl 7ty. Note that because of the different scaling weights carried 
by the frame fields, the components of T A B obtained by functionally differentiating 
A^ A ) do not all have the same eigenvalue A: 

^(A) _^ jnO (A) jnO (A+l-z) rpl (A+z-l) jnl (A) /gg\ 

Similarly, the fact that if) carries a non-trivial scaling weight implies 

A^vr^. (60) 

Expanding (157|) in dilatation eigenvalues, we have 

(z + 2 - A)A( A ) = zT° {A) + TY A) - A_^ A - A - } (61) 

= 2^° (A) + 2vrV A) + za^ A) + (z - A.)^"^. 

Note that (16T1) does not determine A^ A -* for A = z + 2. This is the leading term 
in the dilatation expansion associated with the response functions; it is also the term 
which makes a finite contribution to the on-shell action, as it scales like r~( 2+2 ) at 
leading order, cancelling the scaling of the integration measure \/—h. From the field 
theory point of view, the expectation values of the operators Tab and O^ are the 
functional derivatives of the on-shell bulk action with respect to the boundary data. 
Since the boundary data give the leading terms in the bulk fields, the functional 
derivatives with respect to boundary data can be exchanged for the ones involving 



7 Note that, as stressed recently in [23| . this statement depends crucially on satisfying our asymp- 
totic boundary conditions, and it is difficult to extend this dilatation approach to consider even 
perturbative violations of the boundary conditions. 

21 



&> A \ ip m ( 152|) . (1531) up to some overall powers of r. Once we have cancelled the 
divergent contributions, the leading contribution to the action comes from A^ +2 \ so 
the expectation values will be given by its functional derivatives, that is, 

(£) = r ^T° iz+2 \ (£<> = r 2*+i r ' o (2*+i) 5 {Vi) = r 3 T o^(3) 5 (n ^ = rZ+2T y z+ 2) ; 

(62) 
and ((Dip) = r A *7Q * (in all these quantities, the limit as r — > oo should be taken to 



compute the actual vevs). Note also that for A = z + 2 the right-hand side of (161D 
must vanish. This is precisely the tracelessness condition for the stress tensor complex, 
generalised to the case where we explicitly break scale invariance by introducing a 
source i/j for O^,. 

Our interest is not however primarily in this undetermined part, but in the terms 
in the dilatation expansion of the on-shell action which are locally determined by the 
boundary data. These are determined by (IBTj) together with the constraint equation 
(134|) . whose expansion in dilatation eigenvalues is 

J2 \2K^K^ - 2K A %K AB ^~^ - tt^V^ - -L^ A n A f\V A 7i A y A -^ 

lit 

(63) 

1 _-CA /iln ../A/91-.lDM/n^ 1 (A/V\ /t/A/nl 1 /KfA/o\o 



s<A/2 



l^(A/2)2 _ K (A/2) K ^(A/ 2 ) _ 1 (A/ 2 ) 7r A(A/2) _ J_ (^4^/2)2 



src- '. 



where src^ is the source contribution from the right-hand side of f l34|) . which only 
appears for some particular values of A. 

We want to manipulate this expression to obtain a formula for the right-hand side 
of (IBTl) . To do so, it is useful to note that the first two terms in (l63l) can be rewritten 
as 2K^K^~^ - 2K A s) B K ABi ~ A - s ^ = -2K ( a %ti ab ^- s \ Now from the above large r 
expansion, we can see that 

K° (fS > = z, KJP^S'j. (64) 

For the vector momentum, we can write 

n A = A A + A B K B A -d A A r , (65) 

which gives 

7Tq = aK q = za (66) 

and a constraint on 717; , 



7T, 



(A_) _ J, , _^0 (A_) , j^O (0) _ _^0 (A_) 



i(j + «F l j + ^ o = «AT o + (- 21 - A -)^- ( 67 ) 



Note that we can combine these results to obtain 

-^ (°) _ o^ i iuA 



r B w = -2(z+i)^ B . (6^ 
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For a generic value of A, the terms in ( 163]) with s = and s = A_ are then 



-2Z7T, 



(A) 



2n 



/ (A) 



-(zT° Q (A) 






,1 (Ah 



(A-) 



+ 2r; ^)-(7r^-^ 



7T 



0(A-2A_) 



Putting this together with ( 1611 gives 



(69) 



( 2 + 2-A)A (A) = zT° ( A) + T/ 



,/ (A) 



A_V< 



(A-A_) 



—quad" ' — 



src 



(A) 



(70) 



where quad^ is the remaining quadratic terms, 



qua 



d^ 



0<s<A/2;s^A_ L 

+ [k {A b ] T A ^ A - A ^ + K {A ~ ] ^ A - 2A -^ + Ttf - ) 7T J ( A - A - 



+ 



^(A/^^BCA/S) + ^ 7r (A/2)^(A/2) + _£_ ( y^^A/^ 



l AB 



2m 2 



(71) 



and the source terms will be given below. This is the equation we wish to solve to 
determine the divergent parts of the on-shell action. The advantage of this dilatation 
expansion, relative to the expansion in powers of r we considered in the previous 
section, is that the curvature terms in ( 134"|) . which involved an infinite series of terms 
in the expansion in powers of r, contribute a finite set of terms with definite dilatation 
eigenvalues, so we can write them all down explicitly, and the sub-leading terms 
contribute to further sub-leading terms only through the quadratic terms ( 17T1) . 

Let us first consider the source terms. For the Ricci rotation coefficients Q AB C , 
using the behaviour of the frame fields, the dilatation dimensions are A = 1 for Q QI ° 
and &ij K , A = z for £l 0I J , and A = 2 — z for Qjj°. For the vector field, ifi has by 
definition dimension A = A_. Hence F Q j = — diA + 2Q 0I °A will contain components 
of dimension 1, 1 + A_, and Fjj = 2Q IJ °A will contain components of dimension 
2 — z, 2 — z + A_. The Ricci scalar in terms of the Ricci rotation coefficients is 



R 



-Ad A n CA n + 4Q CA M B 



c 



C" AB 



O nABC 



(72) 



This will have contributions at dimensions 2, 2z and 4 — 2z. 

Putting all this together, we can obtain the source terms at each order: 



A = 0: 



A = 2: 



src 



(A) = _2A + -mV = 2(z + l)(z + 2). 



(73) 



src (A) = R mJ_ F w F AB(D = - Adl n CI c +4n CI c n B IB -2n 0I0 n 0I0 -n IJK n IJK +n 0I0 n 0I0 a 

(74) 
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A = 2z: 

S rc (A) = r {2z) = -Ad n I0 j + An I0 jn J0 j - 2n I0J n I0J . (75) 

A = 4 - 2z: 

src^ = R^~ 2 ^ - \ F ^) F ^{2~ Z ) = _ QlJoQ uo + niJ0 rfJ0 a 2 (?6) 

A = A_: 

src- ' = m 2 atp (77) 

A = A_ + 2: 

W A ) = ~F^ B F AB{1+A ^ = -(djif; + 2Q 0I °i/;)Q I0 a. (78) 

A = A_+4-2z: 

src (A) = l F (fc*) F AB(*-*+*-) = 2n IJ0 n IJ0 atjj. (79) 

A = 2A_: 



src (A) = _ m 2^2 ^ 

A = 2A_ +2: 
src^ = _l F a+A_) FJ 4 B (i + A_) = _1 ( ^ + 2^)^ + 2Q 0J ^)5 IJ . (81) 

A = 2A_ +4-2^: 



src ( A ) = A f (? b *+ a -) f ab(2- z+ a-) = n IJO tfjQ^ (82) 

We will have contributions from the src^ term in (I7D|) only for these values of A. 

There is one mixed term in (i7T|) which involves derivatives of the momenta, namely 
Vati" A - This can be written as 

\/ A 7lA = d 7T - 2O C0 C 7r + d 7 7i7 - 2fi c/ C 7T J , (83) 

so a term in 7r° with dilatation eigenvalue A will contribute a term in V a ka with 
eigenvalue A + z, while a term in 717 with eigenvalue A will contribute a term in Vatt a 
with eigenvalue A + 1. 

We casn now verify that we can solve f!70|) for A^ A -* at each order in terms of the 
sources and the functional derivatives of A^ A ' at earlier orders A' < A. That is, we 
want to see that no term coming from A' > A can appear in the quadratic term 
( 17T|) at order A. This is not immediately obvious because as noted in (j59|) and (16TJJ) . 
the functional derivatives can give terms with lower A, as X^ determines T T 
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and tt\ . The potential problem with ir^ is easily dealt with; in the explicit 

quadratic term, the sum only involves tTq f° r s > A_, so the contribution to the 
quadratic at order A involves n^ ~ for A' < A. n^ also contributes to V a %a, but 
this contribution involves A' + (z — A„), so since A_ < z for all z, this involves only 
A' < A. For the T 0I piece, when T^ +1 ~ z appears in the quadratic with K^ , this 
will give A' < A if s < z — 1. For z < 2, the smallest eigenvalue in Kqj (coming 
from \( 2 >) is s = 3 — z > z — 1 as z < 2. For z > 2, setting to zero the source terms 
for Si implies (as we will see below) that the smallest eigenvalue in Kqj is s — z + 1. 
Thus, all terms in the quadratic piece in (1TT1) are determined by A^ A ' for A < A', 
and we can solve (ITUl) to determine A^ A ) at each order in terms of the pieces at earlier 
orders. In particular, this implies that none of the X^ for A < z + 2 depend on the 
undetermined piece \( z+2 *> which is a non-local function of the boundary data. That 
is, all the divergences in the naive action (ITUj) are local functions of the boundary 
values of the bulk fields, determined by solving (ITUj) recursively. We can then cancel 
the divergences A" for A < z+2 by subtracting these local functions of the boundary 
values of the bulk fields as a local counterterm in our definition of the action. 

The solution for the on-shell action in this dilatation expansion for arbitrary values 
of the source terms will involve dilatation eigenvalues which are the same as the 
powers of r appearing in the discussion of the solution of the equations of motion in 
the previous section, namely 

A = 2n + 2mz + p(4 - 2z) + qA^, n,m,p,q G N. (84) 

We have p = q = if we set the sources for O^ and £ l to zero. There will be additional 
dilatation eigenvalues which appear in the expansion of the on-shell action if we take 
a non-zero value for the response function at A = z + 2, but the terms which are 
determined in terms of the boundary data will have the eigenvalues 



6.1 Counterterms for z < 2 

Let us now demonstrate this approach by constructing the divergent pieces of the 
action for z < 2, thus identifying the required counterterms. One would like to give 
a comprehensive survey of the divergent contributions; unfortunately however, for 
some values of z, there will be a large number of values of p and q in (|84|) which give 
A < z + 2, so a comprehensive discussion would be long and not very illuminating. 
We will therefore confine ourselves to the terms with p = and q = 0,1. 

Consider first the terms with m = n = 0, which don't involve derivatives of the 
boundary data. The first contribution, at A = 0, is determined by the leading falloff 
of the bulk fields. We have the T B given by fl68|) . which gives the leading term in 
the on-shell action to be 

T (0) rpl (0) 

A(0) = z + 2 =- 2 (* + 1 )- ( 85 ) 

Note that this is consistent with re-obtaining T° and T T j as the functional derivatives 
of A using ( 1521) . These values are also consistent with the constraint (1541) . using the 
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value (175]) for the sources. At A = A_, applying (1701) gives 

A (A " } = -zai/>. (86) 

The functional differential equations (1521) . (|53|) can be checked to be consistently 
satisfied, with T B = — zai/)S B . These first two terms reproduce the no-derivatives 
counterterm given in [21] to the relevant order: there we had 

Sa = 4 + zay/-A a A a = 2(z + 1) + zaip + 0(ip 2 ) . (87) 



It is only the terms at linear order that mattered in [2T], as there we only considered 
divergences in the absence of sources. However, the previous calculation has to be 
corrected at higher orders. For example, at A = 2A_, after a little more calculation, 
we find 

A (2A_ ) = z _( 2z -i-AU 2 . (88) 

z + 1 



This is consistent with the constraint on tcq in ( 1671 ). Note that as the frame 



„(A_) 



fields cannot appear undifferentiated in A, possible terms in T J and T° 7 arising from 
functional derivatives of the no-derivative terms will vanish. 

Let us now consider the derivative terms. For A^ there are no quadratic contri- 
butions, and putting in the source contributions in (I70|) gives 

- z \w = rW~f { 2 b f ab u = -4d I n CI c +4n CI c n BI B -2n 0I0 n 0I0 -n IJK n IJK +2n 0I0 n 0I0 a 2 . 

(89) 
For \( 2z \ there is also a contribution from the quadratic term involving V" 4 ^, and 
we have 

-(2-z)\M = rt 2z) + -^(V A n A p 2 = -49 fi /0 / +4(fi /0 / ) 2 -2fi /0J fi /0J + a 2 (n /0 / ) 2 . 

2m 2 

(90) 

On the intermediate steps, the superscripts (2), (2z) etc indicate that we keep only 
the term of the stated dilatation eigenvalue in the expression. The interesting point 
to note here is that these expressions are all by construction covariant with respect 
to coordinate transformations, but because of the different coefficients on the left- 
hand side, the sum of these two-derivative terms is not covariant with respect to the 
Lorentz transformations acting on the frame indices. This is as we would expect, 
given the different scaling dimensions of the frame fields. Put another way, the bulk 
matter field has singled out a particular direction, so the and I indices are treated 
differently. Here we see explicitly how this enters into the required counterterms in 
the on-shell action. For z = 1, the two terms will combine into a Lorentz-covariant 
expression, and this reduces to the usual AdS counterterm. 

By taking functional derivatives of these terms, we can evaluate the corresponding 
components of the stress tensor. For general X^ A \ the functional variation will be 
complicated to calculate, but because these low-order terms can be written as the 
components of a given dilatation eigenvalue in a boundary scalar, we can simply 
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calculate the variation of the scalar quantity and then extract the term of the correct 
dilatation eigenvalue. For X^, this gives 

zT$ = -[2i2oo + R + F 0C F C + ^F 2 - aV c F%f\ (91) 

zTf) = [-2Rjj + 5jjR - F IC Fj C + ^5 u F 2 f\ 

rp(3-z) _ _op(3-z) 
^0/ ~~ zrL 07 ' 

zT^^i^Roj + aVcFT^- 
While for \( 2z \ we have 

(2 - z)Tg z) = -[2Ro + R + ^Vo(V c A c )p\ (92) 

(2-z)Tfj Z) = [-2R IJ + 5 IJ Rf z \ 

,(24 
01 1 

2, 



(2-z)T^ +1) = -2R^ n 



(2 - z)T%* -V = [-2R 0I - Z -^W ^ A c )f^ l) . 

Derivative terms with q 7^ are more involved; to illustrate this, we consider 
the term at A = A_ + 2, which will also be useful for determining 7A , which will 
contribute to quadratic terms in the next subsection. This is the first term for which 
a quadratic term in the stress tensor appears. We have 

(z - A_)A( A - +2 ) = -Kf B ] T AB ^ + -F^ B F AB ^ +A -\ (93) 

Using the values of T AB ~ and 7Tq obtained above to evaluate K A % , we can rewrite 
this as 

n( A -+ 2 ) = 

2{z + l 

(94) 
This gives 

U-A Jtl ( , 2) = _, a ,, \z(3z - A_)T°°( 2 ) + z(2z - 1 - A_)T 7 / 2) + (z + l)n IJO n IJO 



(z-A^X^-^ = , aij . [2(32 - A_)T°°( 2 ) + z(2z - 1 - A_)T 7 / 2) + (z + lRjo^" 



'V> 



2U+1 



(95) 
We see that already at this two derivative level, the form of the counterterms is 
reasonably complicated. Nonetheless, the conceptual picture is simple: As in the 
AdS case, the divergent terms are written explicitly in terms of bulk fields evaluated 
on the boundary. Thus, adding these same expressions to the naive action f jTQj) will 
cancel the divergences. 

6.2 Counterterms for z > 2 

In the asymptotic analysis, we saw that our spacetime would be asymptotically locally 
Lifshitz for z > 2 only if we set to zero the sources for the matter operator O^ and the 
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energy flux £ l . We will therefore consider here only the counterterms for divergences 
which arise with those sources set to zero. To consider correlation functions involving 
Of or £ l , we would need to consider the corresponding sources at the perturbative 
level, following the interesting recent analysis in the AdS case in [231 El] • We leave 
such consideration for future work. 

Setting these sources to zero implies that the dilatation expansion of the action 
involves only terms with A = In + 2mz. We want to consider the divergences in the 
expectation values, and not just the on-shell action, so we should note that setting 
the source for £ l to zero implies that the terms in T Q j with A = 2n + 1 — z will 
necessarily vanish, as there is no local function of the sources with this dilatation 
dimension when fljj = 0. Thus, in the expansion of Too, Tjj there will only be terms 
with A = In + 2mz, while in the expansion of Th, Tjo there will only be terms with 
A = 2n + l + (2m + l)z. 

Setting the sources for O^ and £ % to zero simplifies the expansion, and we can 
now discuss all the divergent terms: for the action, the divergent terms have A = 2n. 
We obtained A^ -* and X^ in the discussion above; we can now consider for example 
A*- 4 -*, which is the only additional term required for z < 4. There is no source term for 
A = 4, so the relevant term in the on-shell action is just given by quadratic terms, 

{z _ 2 )A( 4 > = - Kf B n AB ^ - ^V^. (96) 

Using the expressions for Tjjj and ti q in the previous subsection, one could compute 
this counterterm explicitly, but the full expression is unenlightening. What's impor- 
tant about this calculation is not to obtain the explicit forms of the counterterms, 
but to see that the divergences can be cancelled by adding appropriate counterterms 
which are local functions of the boundary values of the bulk fields. 

7 Discussion 

In this paper, a geometrically natural definition of asymptotically locally Lifshitz 
spacetimes was given. This gives boundary data which act as sources for the stress 
tensor complex of the dual field theory, and can be interpreted as giving a natural 
non-relativistic curved geometry on the boundary. We have shown that solutions of 
the equations of motion for arbitrary values of the boundary data exist for z < 2 in 
the context of the massive vector theory. For 2 < z < 4, bulk solutions satisfying 
these boundary conditions exist if we set some of the sources to zero. This is natural 
from the field theory point of view, as the sources we are setting to zero correspond 
to irrelevant operators. For z > 4, however, we would have to set one of the response 
functions to zero to satisfy these boundary conditions, and we argued in section 14.11 
that these boundary conditions will have to be appropriately modified. 

For arbitrary boundary data, there are divergences in the bulk action. In section 
|6l we showed that for asymptotically locally Lifshitz boundary conditions, all the 
divergences are local functions of the boundary values of bulk fields. Thus the diver- 
gences can be cancelled by adding appropriate local counterterms to the action. We 
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gave a procedure for explicitly deriving these counterterms, and calculated the first 
few terms explicitly. The boundary conditions and the construction of counterterms 
here have a very similar spirit to the familiar AdS calculations; the main message 
of this work is that dealing with Lifshitz only requires a modest extension of the 
holographic toolbox. 

We should revisit an issue to do with divergences from the previous linearized anal- 
ysis. In [2T], there were divergences in £ l coming from the other response functions 
for z > 2, and derivative counterterms were introduced to cancel these divergences. 
In the general analysis we have undertaken here, however, we found that the only 
divergences in the action came from the terms A^ for A < z + 2, and the required 
counterterm action was uniquely determined as a function of the boundary values 
of the bulk fields by cancelling these divergences. This is not a contradiction; as 
observed in |23J , the counterterm action is a function of the boundary values of the 
bulk fields, not of the boundary data. Evaluated on a cutoff surface of finite r and 
written in terms of the linearised mode solutions, this will have a subleading part 
involving the response functions. When we consider irrelevant operators, [23] showed 
that some of these subleading terms can be divergent. That is, when we write the 
divergent part of the action in terms of the linearised modes, it can appear to have 
a divergence involving the response functions, but when we rewrite it in terms of the 
boundary values of bulk fields, it depends locally on these fields. These were called 
pseudo- non-local divergences in [23] . Note that the reason that the calculation is sim- 
pler here than in [23] is that we have switched off the sources for irrelevant operators, 
so that we can still apply the dilatation expansion technique. 

In the approach taken here, the counterterms are necessarily local functions of the 
boundary values of bulk fields, so they should automatically cancel the divergences 
seen in [21] . This can be verified at the two-derivative level by checking that the coun- 
terterm needed to cancel A^ is essentially the same as the two derivative counterterm 
introduced in appendix A of [21Jj, for the parameter o\ = —-. The calculation there 
then shows that this counterterm will cancel the k 4 divergence considered there. 

The analysis given in this paper is essentially complete for z < 2. For z > 2, our 
restriction to asymptotically locally Lifshitz spacetimes forces us to set sources for 
irrelevant operators to zero. It may be interesting to analyse the sources for these 
irrelevant operators perturbatively, following [23J; this would be important particu- 
larly if there were problems where we wanted to consider finite density for the energy 
flux S l . There is also an issue with the calculation of the expectation value of the 
scalar operator for z > 2; the operator becomes irrelevant in this range, but more 
significantly, the source for this operator appears linearly in the frame fields as well as 
in the vector field. It would be interesting to see if similar issues arise for the matter 
content in models which can be embedded in string theory. If so, it would be valuable 
to understand the correct prescription for calculating this expectation value. 

For z > 4, there is a more substantial issue, as asymptotically locally Lifshitz 
boundary conditions require that one of the response functions vanishes, leaving us 
without sufficient freedom to satisfy boundary conditions in the interior of the space- 



i After correcting a sign error there. 
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time. In this regime, our asymptotically locally Lifshitz boundary condition will need 
to be replaced by a different condition. In section I4.1[ we saw that at the linearised 
level we could either adopt a mixed boundary condition for the frame component e\ , 
which would fix the source for Vi as before, or we could continue to impose a simple 
Dirichlet boundary condition for e t , which would now correspond to fixing Vi, and 
regarding it as a source for a vector operator of dimension z — 1. Understanding the 
appropriate boundary condition at the non-linear level and extending the analysis 
performed here to such boundary conditions is perhaps the most important direction 
for future work here. 
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